arXiv: gr-qc/0112017v3 15 Feb 2002 


Draft version February 7, 2008 

Preprint typeset using lATpX style emulateapj v. 04/03/99 


FREQUENCIES OF F MODES IN DIFFERENTIALLY ROTATING RELATIVISTIC STARS AND 

SECULAR STABILITY LIMITS 

Shin’ichirou Yoshida, Luciano Rezzolla 1 
SISSA, International School for Advanced Studies, Via Beirut 2-4, 34014 Trieste, Italy 

AND 

Shigeyuki Karino, Yoshiharu Eriguchi 

Department of Earth Science and Astronomy, Graduate School of Arts and Sciences, University of Tokyo, Komaba, Tokyo 

153-8902, lapan 

Draft version February 7, 2008 

ABSTRACT 

We have computed the eigenfrequencies of / modes for a constant-rest-mass sequences of rapidly rotating 
relativistic inviscid stars in differential rotation. The frequencies have been calculated neglecting the metric per¬ 
turbations (the relativistic Cowling approximation) and expressed as a function of the ratio between the rotational 
kinetic energy and the absolute value of the gravitational energy of the stellar model f3 = T/\W\. The zeros 
and the end-points of these sequences mark respectively the onset of the secular instability driven by gravitational 
radiation-reaction and the maximum value of /3 at which an equilibrium model exists. In differentially rotating 
stars the secular stability limits appear at a (3 larger than those found for uniformly rotating stars. Differential 
rotation, on the other hand, also allows for the existence of equilibrium models at values of (3 larger than those 
for uniformly rotating stars, moving the end-point of the sequences to larger (3. As a result, for some degrees 
of differential rotation, the onset of the secular instability for / modes is generally favoured by the presence of 
differential rotation. 

Subject headings: stars notation — stars:oscillation — stars:neutron 


1. INTRODUCTION 

Instabilities of rotating relativistic stars have been studied for 
more than three decades in relation to the emission of gravita¬ 
tional waves driven by radiation-reaction. Since the discovery 
of the so called Chandrasekhar-Friedman-Schutz (CFS) insta¬ 
bility (Chandrasekhar 1970 ; Friedman and Schutz 1978), in 
fact, the stability properties of a number of non-axisymmetric 
stellar oscillations have been investigated using different tech¬ 
niques (see Friedman 1998 ; Andersson and Kokkotas 2001 for 
recent reviews). The / (fundamental) modes were among the 
first modes of oscillations ever to be studied in rotating relativis¬ 
tic stars. These are spheroidal modes with harmonic indices 
l = m and represent the generalization of the Kelvin modes 
of Newtonian Maclaurin spheroids to compressible fluid stars. 
Over the years, the literature on the subject has been contin¬ 
uously updated and the limits for the onset of the instability 
have been improved through different approaches. This has 
been done by several groups focussing on Newtonian stellar 
models (Bardeen et al. 1977 ; Comins 1979a, 1979b ; Clement 
1979 ; Durisen and Imamura 1981 ; Imamura et al. 1985 Man- 
agan 1985 ; Ipser and Lindblom 1990) on post-Newtonian stel¬ 
lar models (Cutler and Lindblom 1992) and on rapidly rotating 
relativistic stars (Yoshida and Eriguchi 1997 ; Stergioulas and 
Friedman 1998 ; Morsink et al. 1999 ; Yoshida and Eriguchi 
1999). Part of the interest in these modes comes from the fact 
that they have long been regarded as the modes of oscillation 
most susceptible to the CFS instability. In particular, for stars 
rotating at a rate at which mass starts being shed at the equa¬ 
tor (the mass-shedding limit), general relativistic calculations 
on the stability limit indicate that the l = m = 2, /-mode 
(also referred to as the “bar-mode”) could have the shortest 
growth timescale. Because of this, and because of a possible 


weakening of the bulk viscosity at very high temperatures (Lai 
2001), the bar-mode may represent the most important non- 
axisymmetric instability in very hot and rapidly rotating newly 
born neutron stars. 

The occurrence of an /-mode instability in a newly born 
uniformly rotating neutron star is in general prevented by the 
large bulk and shear viscosities of nuclear matter in those con¬ 
ditions and detailed calculations have shown that the instability 
is suppressed except for very large rotation rates, close to the 
mass-shedding limit (Ipser and Lindblom 1991); for superfluid 
stars an even stronger damping was calculated (Lindblom and 
Mendel 1995). More recently, however, a number of new ele¬ 
ments have improved our understanding of the instability and 
have again increased the expectations that the / mode instabil¬ 
ity might characterize the earliest life stages of a newly formed 
neutron star. The first of these new elements was provided by 
Cutler and Lindblom (1992) who have shown that, within a 
post-Newtonian approximation, general relativistic effects tend 
to further destabilize the / mode, lowering the critical value 
of the ratio between the stellar rotational kinetic energy and 
the absolute value of the gravitational energy, (3 C = {T/\W\) C 
at which the secular /-mode instability is triggered. Analo¬ 
gous results have also been found by Yoshida and Eriguchi 
(1997) within the relativistic Cowling approximation and by 
Stergioulas and Friedman (1998) in fully general relativistic 
calculations. The second new element was provided by Shibata 
and Uryu (2000) whose fully general relativistic hydrodynami- 
cal simulations have shown that the remnants of binary neutron 
star mergers could be, at least for polytropic equations of state, 
rapidly and differentially rotating stars. In addition to this, Liu 
and Lindblom (2001) have recently computed the structure of 
objects formed in accretion-induced collapse of rotating white 
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dwarfs and found that these objects can rotate extremely rapidly 
and differentially. Furthermore, a differentially rotating object 
could be produced also during the process of core-collapse in 
supernova leading to the formation of a neutron star. 

Within this framework, differential rotation has two major 
consequences. Firstly, it allows for the existence of an equilib¬ 
rium model at values of f3 which are considerably larger than 
the ones supported by the counterparts with the same rest-mass 
but uniform rotation (Baumgarte et al. 2000). Secondly, as 
shown by a number of authors for Newtonian stellar models 
(Managan 1985, Imamura et al. 1985, Imamura and Durisen 
2001), it increases the critical value (3 C for the onset of the sec¬ 
ular instability. In this revised picture, we have computed the 
eigenfrequencies of / modes and determined the secular stabil¬ 
ity limits for rapidly rotating relativistic stars with differential 
rotation. 

2. FORMULATION 

The equilibrium stellar models are assumed to be stationary 
and axisymmetric and are constructed with a numerical code 
based on the method of Komatsu et al. (1989a,b). Their space- 
time is therefore described by the line element 


Here Ho is the angular velocity at the centre of the coordinate 
system, i.e. at the intersection of the rotation axis with the equa¬ 
torial plane of the equilibrium models, and A R is a dimension¬ 
less parameter accounting for the degree of differential rotation. 
In particular, the degree of differential rotation increases with 
A” 1 , and in the limit of A“ 1 —> 0, the profile reduces to that of 
uniform rotation. In the Newtonian limit, the differential rota¬ 
tion law (jlj) reduces to the so-called “j-constant” law (Clement 
1978) 

- = Al . (4) 

Ho A 2 + r 2 sin 2 9 

and is commonly used in Newtonian calculations (Eriguchi and 
Muller 1985 ; Karino et al. 2001 ; Rezzolla and Yoshida 2001). 

Once the rapidly and differentially rotating equilibrium stel¬ 
lar model has been constructed, we introduce “adiabatic” per¬ 
turbations in the fluid variables so that the adiabatic index of 
the perturbed matter coincides with the polytropic exponent of 
the equilibrium configuration 


ds 2 = -e 2u dt 2 +e 2a (dr 2 +r 2 d0 2 )+e 2 ^r 2 sin 2 9(d(l)-u;dt) 2 , 

( 1 ) 

where u, a, //. and u> are the “gravitational potentials” and are 
functions of the r and 9 coordinates only. The stars are modeled 
as relativistic polytropes with the equation of state (EOS) 

e = p + Np , (2) 

where p, p. e are the pressure, the rest-mass density and the to¬ 
tal energy density, respectively. The subscript “c” refers to the 
maximum value for the equilibrium model which is used for 
normalization. In order to investigate how the limits for the 
secular instability depend on the “stiffness” of the EOS, we 
have performed calculations for two different polytropic mod¬ 
els whose properties are summarized in Table I (We adopt units 
in which c = G = M 0 = 1.). Along each sequence of rapidly 
rotating stellar models, the polytropic constant k, the polytropic 
index N and the total rest-mass Mq are kept constant. 


i-<i) 


pV +1/jv 



N 

K 

Mq/Mq 

Model (a) 

0.5 

6.02 x 10 4 

1.60 

Model \b) 

1.0 

1.00 x 10 2 

1.52 


TABLE I: Properties of the polytropic equilibrium models. 

The rapidly rotating stellar models are constructed after spec¬ 
ifying a law of differential rotation H = H(r, 9) with a choice 
which is, to some extent, arbitrary. This is because there are 
a number of differential rotation laws that satisfy the integra- 
bility condition of the equation of hydrostatic equilibrium as 
well as the Rayleigh criterion for dynamical stability (Tassoul 
1978). All of these laws are physically consistent, cannot be 
excluded on physical grounds and might influence the qualita¬ 
tive behaviour of the results. We here follow the formulation 
suggested by Komatsu et al. (1989a,b) who have modeled the 
rotational angular velocity profile as 


A 2 (fi 0 -fi) = 


(H - oj)r 2 sin 2 9e 2 ^~ v - 
1 — (H — u>)r 2 sin 2 9e 2 ^~ 


(3) 


A p 
P 



Ap 

P 


(5) 


Here, A refers to a Lagrangian perturbation with displacement 
vector £ and is related to the Eulerian perturbation 6 of a generic 
quantity $ through the relation: AT = (-)<!> T £|4>, with £ - be¬ 
ing the Lie derivative along £ (Friedman and Schutz 1978a). 
Because of the stationary and axisymmetric background, the 
Eulerian perturbations are naturally decomposed into a har¬ 
monic component of the type ~ exp (—iat + im(j>), where er 
is the mode angular frequency and m is an integer. Working 
within the Cowling approximation (McDermott et al. 1983), 
we do not include Eulerian perturbations in the metric and the 
fluid oscillations are investigated in the fixed background space- 
time of the equilibrium star (jlj) 2 . This simplification clearly 
introduces an error which is particularly large in the case of 
the lower order / modes. However, as discussed in previ¬ 
ous works (Lindblom and Splinter 1990 ; Yoshida and Ko- 
jima 1990 ; Yoshida and Eriguchi 1997, 1999), the results ob¬ 
tained with this approximation reproduce well the qualitative 
behaviour of the mode eigenfrequencies with the errors being of 
the order of 10% for low mode-numbers and progressively less 
for high mode-numbers (Lindblom and Splinter 1990 ; Yoshida 
and Eriguchi 1997, 1999). Furthermore, by neglecting the con¬ 
tributions coming from the metric perturbations, the Cowling 
approximation tends, at least for the lowest mode-numbers (i.e. 
to = 2,3), to overestimate the stability (see the comparison in 
Figs. 2 and 3). 

Introducing the harmonic perturbations in the hydrodynam- 
ical equations yields a system of four partial differential equa¬ 
tions accounting for baryon number conservation and conserva¬ 
tion of the stress-energy tensor. These equations are then solved 
for the four unknown functions represented by the components 
of the Eulerian 3-velocity perturbation 5v r , Sv e , Sv^ and the di¬ 
mensionless quantity q = dp/(t + p). The partial differential 
equations are discretized on a two-dimensional numerical grid 
and are solved following the strategy discussed by Yoshida and 
Eriguchi (1999), after imposing Ap = 0 at the surface of the 
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unperturbed star as a boundary condition. For each degree of 
differential rotation .4 n and for each mode-number to, the solu¬ 
tion to the eigenvalue problem is found for increasing values of 
/3 until the limit of mass-shedding is found, which is indicated 
with (3 S (see Komatsu et al. 1989a for a definition of (3). 

3. RESULTS 

The results of these calculations are presented in Fig.fi where 
we plot the frequencies of the m = 2 mode as a function of the 
parameter 3 for the model (6) of Table I. Different curves re¬ 
fer to different degrees of differential rotation. Note that the 
values of (3 at which er = 0 (i.e. (3 C ) signal the onset of the 
secular instability (neutral points) and that these increase as the 
degree of differential rotation is increased. At the same time, 
however, the differentially rotating models are able to support 
larger values of (3 before reaching the mass shedding limit (3 S , 
represented by the end-points of the curves in Fig ^3j 



Fig. 1.— Eigenfrequencies of the m = 2 mode as a function of the parame¬ 
ter (3 = T/\ W\ for the model ( b ) of Table I. Different curves refer to different 
degrees of differential rotation, with the A" 1 = 0.0 line being the one of a 
uniformly rotating model. The filled dot indicates the neutral stability point of 
a uniformly rotating star computed in full General Relativity (Stergioulas and 
Friedman 1998). 

Given the results of Fig.fi it is natural to ask whether differ¬ 
ential rotation favours or not the onset of the /-mode instabil¬ 
ity in rapidly rotating relativistic stars. The evidence that f3 c 
increases with increasing differential rotation is not sufficient 
to draw a conclusion since differentially rotating models can 
support values of (3 which are considerably larger than those 
of uniformly rotating stars before reaching the mass-shedding 
limit. This is simply due to the fact that in the differentially 
rotating model described by eq. (fi) the inner regions can rotate 
rapidly while the outer regions rotate more slowly, preventing 
mass-shedding. Thus, to quantify the importance of differen¬ 
tial rotation for the onset of the instability we need to measure 
not only the secular stability limit [3 C , but also how close the 
latter is to the ultimate limit of mass-shedding (3 S . A relative 
measure of the two quantities across sequences with different 
degrees of differential rotation can provide information on the 
likelihood of the existence of a configuration at the onset of 
the secular instability and clarify the role played by differen¬ 
tial rotation. In this sense, the normalization of (3 C with 3 S is 


equivalent to the normalization, commonly used for uniformly 
rotating stars, of writing the stellar angular velocity in terms of 
the mass-shedding angular velocity. 

In Fig.fi] we show the behaviour of the ratio (3 C = (3 C /(3 S as 
obtained for degrees of differential rotation ranging between 
A^ 1 = 0 (uniform rotation) to Af 1 = 1.25. To avoid con¬ 
tamination of our results from configurations that start show¬ 
ing a distinct toroidal topology, we stop our sequences for 
A~ x = 1.25 but were able to find eigenfrequencies also for 
larger values of A” 1 (see Table II). 



Fig. 2.— p c as a function of the rate of differential rotation. The two panels 
refer to the two models of Table I and the different line types refer to different 
mode-numbers. The filled symbols show a comparison with the fully general 
relativistic calculations of Stergioulas and Friedman (1998). 

The two panels of Fig.fi] refer to the models (a) and ( b ) of 
Table I, respectively. In each panel, the different lines refer to 
the mode-numbers to = 2, 3 and 4. As anticipated above, the 
Cowling approximation overestimates the stability of the stellar 
models by introducing an error which is of the order of 10% for 
the to = 2 mode but is considerably smaller for higher mode- 
numbers. This can be appreciated by comparing with the fully 
general relativistic values of / 3 C computed for uniformly rapidly 
rotating stellar models (Stergioulas and Friedman 1998) and in¬ 
dicated with filled symbols in the lower panel of Fig.fi The 
overall decrease of (3 C for increasing degrees of differential ro¬ 
tation provides a clear indication that the onset of the secular 
/-mode instability is in general favoured by the presence of 
differential rotation. This is most evident in the case of a stiffer 
EOS [model (a) of Table I] where, for reasonable degrees of 
differential rotation, f3 c is reduced by about 17%. A direct com¬ 
parison with corresponding values in the literature is presented 
in Table II, where we report the values of / 3 C for different mode- 
numbers and degrees of differential rotation. 

To validate our relativistic results, we have computed the 
eigenfrequencies and the secular stability limits for / modes 
in rapidly and differentially rotating Newtonian stellar models 
(A more extended discussion of these results will be presented 
in Karino et al. (2002). The numerical strategy followed is the 
same as the one discussed above for the relativistic Cowling 
approximation but we have considered both the case in which 
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perturbations in the gravitational potential are neglected (New¬ 
tonian Cowling approximation) and the case where they are in¬ 
cluded. The stellar models used in this case are also poly tropes 
with index N = 0.5. The results found for the correspond¬ 
ing Newtonian quantity (/3 C ) N = (/3 C //3 S ) N are presented in 
FigJ^ where we show (/3 C ) N in the Newtonian Cowling approx¬ 
imation (dashed lines) and in the fully Newtonian calculations 
(solid lines), for the m = 3, / mode. As expected, the results 
found in Newtonian gravity are qualitatively similar to the rela¬ 
tivistic ones and show that the presence of differential rotation 
introduces a general decrease of (/3 C ) N . 


V 

TO = 2 

TO = 3 

TO = 4 

0.00 

0.1300 

0.0748 

0.0516 


0.0923 

0.0627 

0.0436 

0.25 

0.1310 

0.0761 

0.0528 


0.0943 

0.0635 

0.0441 

0.50 

0.1410 

0.0798 

0.0552 


0.1060 

0.0668 

0.0471 

1.00 

0.1530 

0.0845 

0.0587 


0.1250 

0.0753 

0.0513 

1.25 

0.1630 

0.0890 

0.0615 


0.1420 

0.0818 

0.0552 

1.43 

0.1710 

0.0927 

0.0641 


0.1580 

0.0865 

0.0583 


TABLE II: /3 C for different degrees of differential rotation. The re¬ 
sults are presented for m = 2,3,4 and the two values for each mode- 
number and A“ 1 refer to models (a) and ( b ) of Table I, respectively. 

There are three aspects of Fig.^| that are worth commenting 
on. Firstly, it should be noted that we have not plotted the be¬ 
haviour for the to = 2 mode. This is because this mode is 
never found unstable within the Newtonian Cowling approxi¬ 
mation and in the range of differential rotation considered here. 
Secondly, when compared with Fig. the difference between 
the two approaches appears to be larger for the Newtonian star. 
This is to be expected since the Cowling approximation is more 
accurate when the longitudinal gravitational field is stronger 
and hence in General Relativity. A similar improvement of 
the Cowling approximation is also observed when the poly¬ 


tropic index is varied. In the case of a softer EOS, in fact, the 
star is generally more compact and the Cowling approximation 
more accurate (Yoshida and Kojima 1997). Thirdly, it should be 
noted that the inclusion of the perturbations in the gravitational 
potential lowers systematically the values of (/? C ) N for the var¬ 
ious degrees of differential rotation indicating that the Cowling 
approximation tends systematically to overestimate the stabil¬ 
ity limits for the onset of the secular /-mode instability. This 
behaviour, encountered also in relativistic gravity, supports our 
expectations that a fully general relativistic analysis, in which 
the metric perturbations are not neglected, will confirm our re¬ 
sults and provide even more promising estimates for the secular 
stability limits of f modes to the CFS instability. 



Fig. 3.— (/3c) n as a function of the Newtonian degree of differential ro¬ 
tation. The different line types refer to the Newtonian Cowling approximation 
(dashed lines) and to full Newtonian calculations (solid lines). 
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